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Abstract. We explore consistency among different distance measures, including 
Supcrnovae Type la data, measurements of the Hubble parameter, and 
determination of the Baryon acoustic oscillation scale. We present new constraints 
on the cosmic transparency combining H(z) data together with the latest 
Supernovae Type la data compilation. This combination, in the context 
of a flat ACDM model, improves current constraints by nearly an order of 
magnitude although the constraints presented here are parametric rather than 
non-parametric. We re-examine the recently reported tension between the Baryon 
acoustic oscillation scale and Supernovae data in light of possible deviations from 
transparency, concluding that the source of the discrepancy may most likely be 
found among systematic effects of the modelling of the low redshift data or a 
simple ~ 2-cr statistical fluke, rather than in exotic physics. Finally, we attempt 
to draw model-independent conclusions about the recent accelerated expansion, 
determining the acceleration redshift to be z aC c = 0.35^q'j3 (1-<t). 



1. Introduction 

In the -extremely successful- standard LCDM model there are well-defined 
relationships between different distance measures such as luminosity distance, 
angular diameter distance, BAO scale, etc. While combining these measurements 
helps breaking parameter degeneracies and constraining cosmological parameters, 
comparing them helps constraining possible deviations from the assumptions 
underlying the standard cosmological model or exotic physics. This is what we set 
out to do here. 

We will first consider cosmic "transparency". The Etherington relation [11] 
provides a distance duality which holds for general metric theories of gravity in 
any background: conservation of phase-space photon density and Lorentz-invariance 
implies that the luminosity distance is (1 + z) 2 times the angular diameter distance, 
where z denotes the redshift. This duality can be a powerful test of exotic physics 
as well as simpler astrophysical effects. For example, it depends crucially on photon 
conservation which can be violated by any source of attenuation. There are several 
known sources of attenuation, which are expected to be clustered e.g., interstellar 
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dust, gas and plasma around galaxies. More speculatively, if dark matter is axion or 
axion-like it can have interactions with photons and produce effective absorption. The 
effect of clustered sources of attenuation, which are then expected to be clustered with 
the large-scale structures, can be constrained with "angular difference" measurements 
and/or with cross correlation with large-scale structures [13 [4] and are at the 0.1% 
level. In principle, there could also be an unclustered component of attenuation, 
which, however, is much more difficult to constrain. For example, there have been 
several attempts to explain the observed dimming of Type la Supernovae without 
resorting to cosmic acceleration, but invoking unclustered sources of attenuation such 
as photon-axion mixing J7] [20] [5] or gray dust [1] . Both of these effects would have 
to have little frequency dependence above microwave frequencies to evade detection 
through reddening while still contributing to overall opacity [22l [23] . 

Radial and angular distance measures can be combined to constrain this 
opacity [2T1 |2l E] ■ In particular, [3T] combine Baryon Acoustic Oscillations (BAO) 
and Supernovae type la data to constrain the difference in opacity between two 
different redshifts using present data and forecast constraints from forthcoming 
BAO measurements. Here we build upon [21]: we combine Supernovae data with 
independent determinations of the cosmic expansion history H(z) [30] to improve 
constraints on possible deviations from transparency. 

Within the ACDM model, this data set combination, even allowing for opacity, 
yields predictions for the observed BAO scale. We next compare these predictions 
with the measured values at z — 0.35 and z — 0.2 10, 25]. [25] has hinted at a 2.4-<r 
discrepancy between distance measures derived from Supernovae and BAO within the 
ACDM model. We revisit here this issue comparing distance measures derived from 
CMB 16 1, from BAO [13 US], and from our own transparency analysis. 

Finally, after exploring the complementarity and consistency of different distance 
measures, we attempt to draw some model-independent conclusions about the recent 
accelerated expansion by combining independent data-sets that are highly consistent 
with each other. In particular, combining expansion history [HI [30], SNIa [24] and 
HST key 113] data, we find evidence for recent acceleration at z < 0.5 and indications 
of past deceleration at z > 1. We determine the transition rcdshift from deceleration 
to acceleration to be z acc = 0.35lg '^3 at the 1-tr level. 

2. Transparency 

Luminosity distance and angular diameter distance are related by the "Etherington 
relation" [TTj : 



This relation holds for general metric theories of gravity in any background; it depends 
only on conservation of phase-space density of photons (transparency) and Lorentz 
invariance. Models that violate Lorentz invariance have non-trivial dispersion relations 
in the vacuum. This effect however generally becomes larger with increasing energy 
and is expected to be immeasurably small at low-energy (visible bands). Transparency, 
however, could be violated by any source of photon attenuation. While angular 
difference measurements [191 [4] have imposed strong constraints on clustered opacity 
(i.e. strong upper limits on differences in opacity along different lines of sight), 
there could still be an unclustered source of attenuation, which is much harder to 
constrain. Examples of effects that appear as uniform attenuation include, gray 



d L (z) = (l + z) 




(1) 
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dust PP, replenishing dust [57] (used to try to explain Supernovae results without 
resorting to cosmic acceleration, or to relax the constraints on cosmic acceleration) 
or more exotic physics such as axion-photon mixing [7J [5]- [2] propose to combine 
Supernovae measurements of with various estimates of dA to constrain exotic 
physics violating Etherington's distance duality. [21] propose to combine BAO and SN 
la data and constrain the difference in opacity between two different redshifts (z = 0.2 
and z = 0.35) from current data. 

In this section we develop this approach further and use H{z) data [HI [30] 
in combination with Supernovae data [17 to constrain possible deviations from 
the "Ethcrington relation" . If there was a source of "photon absorption" affecting 
the universe transparency, the distance modulus derived from Supernovae would be 
systematically affected, in particular any effect that reduce the number of photons 
would dim the Supernovae brightness and increase d^. Let r(z) denote the opacity 
between an observer at z = and a source at z due to e.g. extinction. The flux 
received from the source would be attenuated by a factor exp(— t(z)) and thus for the 
luminosity distanc^l 

dt t obs — d 2 L,true eT ^ ■ (3) 

Therefore, the inferred ("observed") distance modulus would be: 

DM obs (z) = DMtrue(z) + 2.5[loge]r(z) . (4) 

Measurements of DM b s are taken from the Supernovae Union compilation 17J. We 
shall compare this with the -unabsorbed- luminosity distance durue inferred from the 
H(z) measurements of [30) . This measurement is obtained from ages of old passively 
evolving galaxies: it relies on the detailed shape of the galaxy spectra but not on 
the galaxy luminosity. It will therefore not be affected by a non-zero r(z) since r is 
assumed (and constrained by independent observations [2U [53]) n °t to be strongly 
wavelength dependent in the optical band. 

In particular, for a general FLRW cosmology: 

1 .. — T dz' 

^L,true [ 



-) (5) 



HoV^h VV Jo E(z) 
where Sk(x) stands for sin(x), x, or sinh(x) depending on 0^. being positive, zero or 
negative respectively, and 

E(z) = H(z)/H = [a m (l + zf + n v (l + z) 3 ( 1 +»)] 1 /2 . (6) 

To be able to use the full redshift range of the available data, we consider 
the following simple parameterization of a deviation from the Ethcrington relation 
dh = rfyt(l + z)( 2+<i \ with e parameterizing departures from transparency. To 
understand the physical meaning of a constraint on e we note that for small e and 
z < 1 this is equivalent to assuming an optical depth parameterization t(z) — 2ez or 
r = (1 + z) a — 1 with the correspondence a — 2e. While this identification is based 
on a Taylor expansion we have verified that the expansion is good to better than 20% 
for the entire range of e and the redshift range considered and, as explained below, 
to better than 1% in the allowed rage. We consider a flat ACDM underlying model: 

\ If one wanted to interpret the optical depth in terms of the comoving number density of absorbers 
n(z) and their cross section o"(z) it will be 



t(z) = n(z)o-(z)c dz= n(z)o-(z)— ——dz. 

Jo H(z) J H E{z) 



(2) 
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H(z) and thus d^true depends only on f2 m , Hq, while rf_L i0 b s and thus DM b s depends 
on Q m , H and e. We fit the two data sets [3U] and [T7j both separately and jointly, 
imposing a Hubble constant prior from the HST key project [13] • In effect, the HST 
key project introduces one extra data point at low redshift. 

Figure [1] shows the constraints in the O m ,e plane after marginalizatiorj| over Hq. 
The dark blue regions are the 1 and 2-er joint (2-parameters) confidence levels for the 
Supernovae data. The dark blue contours in the background show the constraints 
considering only Supernovae at z < while the dark blue contours in the foreground 
correspond to all redshifts . With the extra degree of freedom of opacity, Supernovae 
data can only impose an upper limit to f2 m even if the universe is assumed to be 
spatially flat. Note that negative values of e would correspond to Supernovae being 
brighter than expected, which is unphysical if interpreted in terms of transparency, 
but it can still be interpreted in terms of departures from the Etherington relation. 
The supernova-only constraints do not extend to arbitrarily large values of e. This is 
due to the presence of high redshift data breaking the degeneracy: the high redhsift 
supernova data show that the Supernovae dimming is evident at low redshift but not 
at high redshift \\7 ) 26\. The lighter blue regions show the corresponding constraints 
from H(z) data. As the H{z) measurement is not affected by transparency, this data 
set yields no constraints on e. The dotted and solid black lines show the combined 
constraints for the case of Supernovae data only at z < 0.5 and for all redshifts 
respectively. 

Figure [2] shows the constraints on e after marginalization over the other 
parameters. We obtain e = — 0.08^0.20 f° r Supernovae data alone (all redshifts) and 
e = —0.01^0 09 when adding H{z) data at 95% confidence. 

If we make the identification of r ~ 2ez, this represents an improvement of almost 
a factor 7 on the constraints obtained in [21] of At = r(0.35) — t(0.2) < 0.13 at 95% 
confidence: we obtain At < 0.02. This can be understood because they used only 
data at two redshifts to constrain At in a parameterization-independent way, while 
we use a range of redshifts to obtain a parameter-dependent constraint on t{z). In 
both analyses one parameter is effectively measured, in the case of [21] the parameter 
is At, while in our case it is e. 

Since our constraints are parametric, we should examine the extend to which our 
results depend on the chosen parameterization. Opacity should increase monotonically 
with redshift, while the explanation of the observed behaviour of supernovae brightness 
as function of redshift in terms of imperfect transparency would require a change of 
monotonicity at z > 1 (supernovae get dimmer at z < 0.3 and brighter at higher 
redshifts). Thus, for generic parameterizations respecting the above monotonicity 
one expects to obtain strong constraints on cosmic transparency, even though the 
quantitative details will depend on the particular parameterization. Further, since 
the H(z) data are not sensitive to opacity, the combined data-set constraints will be 
even less sensitive. 

For example for z < 0.5 we find that e = 0.9 mimics to better than 4% the 
replenishing dust model shown by [57] to explain Supernovae results without resorting 
to cosmic acceleration. Still, Fig. 1 shows that the combined supernovae+i?(z) 
constraints do not degrade significantly. To demonstrate this further we consider an 
alternative parametrization, which is linear in cosmic time rather than redshift, and 

§ This should not be confused with maximization: we do not choose the values of Ho that maximize 
the likelihood but actually integrate the likelihood over the Hq. 
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Figure 1. Confidence levels (68% and 95%) in the O m ,e plane from the joint 
analysis of Supernovae and H(z) data, after marginalization over Ho- Dark blue 
regions are constraints from Supernovae data (contours in the background use only 
data at z < 0.5 while contours in the foreground use all data), light blue regions 
from H(z) data, and black transparent contours are the combined constraints. 
Here we have assumed an underlying flat ACDM model, where Supernovae- 
inferred luminosity distances are affected by a deviation from the Etherington 
relation parameterized by e. 




Figure 2. x 2 f° r e (having marginalised over f2 m and Hq) in the case of SN 
data only (left panel) and using both datasets (right panel). The dashed line 
corresponds to 95% CL, Ax 2 = 4. 



which attempts to model supernova evolution. In this model the apparent Supernovae 
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dimming is given by [12] : 

Am = K^—t- . (7) 

to — t\ 

Here to, t\ is cosmic time at z = 0, 1 respectively, while K is the "opacity" parameter 
to be constrained. In this case we find the constraints to be similar to the dark 
background contours of Fig. 1 : empirically we find that there is a linear correpondance 
between e and K given by e = 0.65K + 0.03 that produces the same contours in Fig 
1. One thus obtains from the joint analysis with H (z) data K < 0.3. 

Further, one may worry that the Taylor expansion we used to connect our 
constraints on e to constraints on t may break down. We have verified that the 
expansion is good to 20% for the entire range of e and the redshift range considered. 
The maximum deviation occurs for e < —0.4 at large redshifts z ~ 1.5, while at 
e ~ —0.08 (constraint obtained from SN data only), and e ~ —0.01 (constraint 
from the combined SN and H(z) data), the only parameter regions where we have 
interpreted e in terms of opacity r, the expansion is accurate to better than 10% and 
1% respectively, even at the maximum redshift considered. 



3. Consistency with the measured BAO scale 

Having allowed the Supernovae data the extra degree of freedom of violating photon 
conservation, we can check the cosmological constraints one can obtain in combination 
with H (z) data. This is shown in figure [3j where the dark blue regions are the 
Supernovae constraints, the light blue ones are for H(z) and the black contours are 
the joint constraints. Note that as seen before the Supernovae data only give an 
upper limit to f2 m and, by construction, supernovae data do not constrain Hq. These 
constraints are in good agreement with the WMAP5-determined parameters [16] . 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3. Likelihood contours for flat KCDM models in the Q m -Ho plane. Dark 
blue regions show the Supernovae constraints (data from all redshifts), while the 
dark blue ones are for H (z) constraints. Black contours show the joint constraints. 
By construction supernovae data do not constrain Hq. 
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A hint of a discrepancy at the 2.4-a level between Supernovae and BAO data at 
redshifts z = 0.2 and z = 0.35 was indicated in [55]. On the other hand, [TUj find 
agreement for BAO data at z = 0.35. To investigate whether this discrepancy can 
be explained by a non-zero e, we set out to compare these measurements of the BAO 
scale with our distance measurements of §1. 

BAO surveys in principle could measure independently H(z) and cIa(z), or, more 
precisely, d,A(z)/r s and H(z)r s from transversal and radial clustering respectively. 
Here r s denotes the sound horizon at radiation drag. Current surveys, however, can 
only measure an angle-averaged distance measure Dy(z) in combination with r s . 



The sound horizon r s is determined exquisitely by CMB data so it is possible 
to translate constraints on r s /Dy into constraints on Dy and vice versa without 
degrading the error-bars. In particular, [10] reports constraints on LV(0.35) obtained 
from SDSS Luminous red galaxies, while |25j reports joint constraints on r s /Dy (0.2) 
and r s /Dv(0.35) from a combination of the SDSS and 2dFGRS surveys. On the 
CMB side, WMAP 5-year data analysis [TB] reports constraints on r s , r s /Dy (0.2), 
r s /Dy (0.35). Finally, our joint analysis of Supernovae and H(z) (Fig[3]) also yield 
joint constraints on Dy (0.2) and IV (0.35). 

In figure S] we show constraints from these data sets separately in the -Dy (0.2), 
-Dy(0.35) plane (top panel) and r s /ZV(0.2), r s /Dv (0.35) (bottom panel). As before, 
the underlying model is assumed to be a flat ACDM. In this figure the transparent 
contours are the constraints obtained from Supernovae and H(z) data allowing for 
non-zero e (corresponding to the solid black lines of Fig. [3]). The thin filled contours 
are the WMAP5-years constraints. The light blue region corresponds to the 1-er 
constraint from [10] (note that [10] constrains only Dy at z = 0.35), and the dark 
blue confidence regions correspond to the 1 and 2-cr joint constraints from [2"5] . 

Figure 2] shows that imperfect transparency cannot explain the tension between 
BAO and Supernovae distances. When allowing for imperfect transparency, our 
joint SN-H(z) constraints give rise to much broader contours than the corresponding 
WMAP constraints, but the reported tension remains even with this different data 
set. Supernovae distances are consistent with both CMB data and high-redshift BAO 
constraints, and only the constraints of 25j show some tension. In interpreting this 
figure one should keep in mind that the data of [25] and [10] are not fully independent: 
the data at z = 0.35 are largely common to the two analyses while the error-bars 
include statistical errors. In agreement with [25j we see that the discrepancy from the 
ACDM model comes mainly from the low redshift BAO distance measure. The value 
of Dy(0.2) appears to be too low by 5 to 10%. Even considering only the Supernovae 
and BAO data of [25] , an imperfect transparency (e > 0) would make the discrepancy 
worse, as Supernovae would favor less acceleration while the low redshift BAO data 
require stronger acceleration. An e < would achieve this, however, it cannot be 
interpreted in terms of reduced transparency, and the addition of H (z) data puts a 
strong lower limit on e, excluding this possibility. It is well known that non linearities, 
(increasingly more important at low redshift), may shift the observed position of the 
BAO as traced by galaxies (e.g., [29] and references therein) and it does so in the 
direction of reducing the inferred Dy. However, this effect is bound to be at, or 
below, the 1% level for the robust technique to recover the signal used by [25]; it 
cannot therefore be solely responsible for the shift in the measured Dy(0.2). Figured] 




(8) 
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Figure 4. Constraints from CMB, BAO, and Supernovae+/i"(z) in the Dy(0.2), 
ZV(0.35) plane (top panel) and r s /D v (0.2), r s /D v (0.35) (bottom panel). The 
underlying model is assumed to be a flat ACDM. The transparent contours are 
the constraints obtained from Supernovae and H(z) data allowing for non-zero e, 
the thin filled contours are the WMAP5-years constraints (1 and 2-cr joint), the 
light blue region corresponds to the 1-cr constraint from |10| (constrains only Dy 
at z = 0.35), and the dark blue confidence regions correspond to the 1 and 2-<r 
joint constraints from 1251 . 



and these considerations seem to indicate that the source of the discrepancy may most 
likely be found among systematic effects of the modelling of the low redshift data or 
a simple ~ 2-cr statistical fluke, rather than in exotic physics. 
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Figure 5. Top: Our data set combination: HST Key project determination 
of Ho (filled square), points from Supernovae I24I (filled dots), points from H(z) 
determinations of I30I ■ Note h(z) = H(z)/100. For comparison two ACDM model 
lines are plotted. Bottom: Quadratic, cubic and quartic polynomial fits. 

Having explored the complementarity and consistency of different distance 
measures and having found that Supernovae and H(z) data are consistent, we now 
attempt to draw some model-independent conclusions about the recent accelerated 
expansion by combining them. There is evidence for recent acceleration at high- 
statistical significance, but most analyses in the literature are done in a model- 
dependent, parametric way, where the acceleration is due to a dark energy component 
with properties specified by few parameters (e.g., f2y, w or fly, w , w a ). A fully non- 
parametric analysis is probably impossible, but we can attempt a model-independent 
approach. 

To do so we follow [HGITJISS] and introduce the function f(z) = (1 + z)/H{z) = 
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a : f(z) has positive slope if the universe is accelerating and negative if it is 
decelerating. In fact the deceleration parameter q = —da/a 2 is related to f(z) by 
q=-H(z)f(z). 

While the H(z) data of [3D] and [TSJ for z = yield f(z) directly, Supernovae- 
determined dh does not, and f(z) can be inferred only taking a derivative [28l I3TI [9] . 
|24] presented such a model-independent analysis and in particular provided estimates 
for f(z) from the same Supernovae data set used here, which we now combine with 
HST-key [13j and determinations of H (z) [30]. The details of obtaining the expansion 
history H(z) by differentiating Supernovae data can be found in Refs. [3TJ[21], the 
method assumes a spatially flat universe so we will only consider the zero-curvature 
case. 

The resulting data points are shown in the top panel of Fig. [5J where two ACDM 
models are also shown for comparison. A visual inspection of Fig. [5] may indicate that 
the data points follow the expected trend (acceleration at low z and deceleration at 
large z), but that given the large errors in /(z), the data may also be consistent with 
a constant function of z (no acceleration/deceleration). 

To investigate the statistical significance of a possible deviation from a constant 
f(z), we consider three different analyses, namely: i) weighted fits with linear, 
quadratic, cubic and quartic polynomials, ii) weighted linear fits up to different upper 
redshifts, and Hi) piecewise linear fits. These should be seen as model-independent 
parameterizations of the behavior of f{z). We use more than one parameterization to 
assess the robustness of our findings: 

Weighted fits with quadratic, cubic and quartic polynomials yield curves that 
are increasing -corresponding to acceleration- up to redshift z acc — 0.75, 0.45, 0.45 
respectively and decreasing -corresponding to deceleration- at z > z acc , although at 
z > 1.5 the reconstruction become noisy. This is shown in the bottom panel of Fig. 
[5] The best fit chisquares are 10.59, 7.22, 7.17 quadratic (3 parameters), cubic (4 
parameters) and quartic (5 parameters) polynomials respectively. The linear fit has 
a best fit chisquare of 13.93. As there are only 14 data points, a simple chisquare 
analysis indicates that a straight line is a good fit and that there is no evidence for 
acceleration or deceleration. Nevertheless, one should consider the possibility that the 
errors may be over-estimated. We thus attempt to treat the linear and polynomial 
fits as nested models and apply information criteria as model selection tool. 

The Akaike information criteria (AIC) for the linear, quadratic, cubic and quartic 
fits are 17.9, 16.59, 15.22 and 17.17 respectively, while the corresponding Bayesian 
information criteria (BIC) are 19.2, 18.5, 17.77 and 20.36. Both AIC and BIC favour 
the cubic fit. In particular, the BIC provides "positive evidence" (in the Jeffreys 
scale |14j ) in favour of the cubic polynomial with respect to the linear model, with a 
Jeffreys factor of A(BIC)^ 2. We interpret this as slight positive evidence for model- 
independent detection of recent cosmic acceleration at redshifts smaller than 0.5 or 
so, and past deceleration at large redshifts {z > 1). It is worth keeping in mind that 
although both relative measures favour the cubic fit, the reduced chisquared for the 
linear model is close to unity suggesting that the linear model is also a good fit to the 
data. We interpret this as indication that the error bars for the used measurements 
may have been overestimated. 

We then try to understand better the result of the linear fit model by considering 
cases ii) and Hi): weighted linear fits up to different upper redshifts, and piecewise 
linear fits. The weighted linear fits as a function of the maximum redshift considered, 
Zmax, are shown in Fig. [6] lines with z max between 0.3 and 1, have a positive slope, 
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and lines with higher z max have flat or slightly negative slopes. Fig. [5] (bottom panel) 
shows the slopes of the lines (with the corresponding 1-er errors) as a function of z max . 
Thus, by fitting a single line from z — to z max we see evidence for acceleration at 





Figure 6. Top: Linear fits up to redshift z max . Lines broadly fall into two 
categories, lines with a positive slope of 0.5 ±0.15 (z m ax < 1), and lines with flat 
slopes for z max > 1.1. Bottom: Slope as a function of z max . Starting from left 
and moving towards higher z max , more points are included in the fitting. 

low z, but no clear indication of acceleration or deceleration when all redshifts up to 
Zmax > 1 are considered. This analysis cannot provide an estimate for the redshift of 
acceleration onset. 

A bilinear fit on the other hand, can give a model-independent estimate of the 
transition redshift z acc . We perform a continuous piecewise linear fit dividing the 
data in two redshift intervals: from redshift to z* and from z* to z = 1.8. In this 
way z* yields an estimate of the acceleration redshift z acc . As before, the reduced 
chisquares are always < 1: any z aC c gives formally an acceptable fit. However, the 
chisquare changes significantly as a function of the parameters. The top panel of 
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figure [7] shows the best fit (acceleration slope df/dz = 0.63 ± 0.18 and deceleration 
slope df/dz — —0.35 ± 0.08 for the best fit z acc = 0.35). The bottom panel shows 
the chisquare as a function of z acc . Note that while goodness of fit tests considers 
the absolute value of \ 2 an d all values of z acc pass the goodness of fit test, parameter 
estimation procedure only considers \ 2 differences from the minimum. The bottom 
panel of Fig. [7] clearly shows that when looking at A% 2 , the data constrain z acc ; in 
particular, z acc = 0.35^o'i3 a t l _cr level. This may indicate that the error-bars may 
be over-estimated or highly correlated. 

2.51 1 1 ' ' ' 1 ' ' 1 i 1 ' 1 ' ' 1 1 ' ' I 
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Figure 7. Top: The best fit continuous piecewise linear fit. Bottom: Chisquare 
as a function of the transition redshift between deceleration and acceleration z aC c- 
Dashed line indicates a A\ 2 = 1. 

The above analysis has been carried out without reference to any specific model 
of the energy content of the universe, nor to any specific theory of gravity. We 
have assumed, however, homogeneity and isotropy at large scales, so that the FLRW 
geometry can be considered a good approximation. In addition, the method used to 
obtain four out of the fourteen data-points used in this section (those coming from 
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differentiating Supernovae data) assumes spatial flatness. 
5. Conclusions 

Combining different distance measurements is a powerful tool to constrain 
cosmological parameters. Comparing them instead helps constrain possible deviations 
from the assumptions underlying the standard cosmological model or exotic 
physics. We have considered Supernovae data [T7J, direct meauserments of cosmic 
expansion [301 H3] and determinations of the B AO scale [25] [TU] . 

We have explored the possibility of uniform deviations from "cosmic 
transparency" through its effects on distance duality, by parameterizing possible 
deviations from the Etherington relation ([T]). Such deviations from the Etherington 
relation might arise from any uniform and "gray" source of attenuation, from 
astrophysical processes such as gray [T] or replenishing dust [27J (used to try to explain 
Supernovae results without resorting to cosmic acceleration, or relaxing the constraints 
on cosmic acceleration), to more exotic physics such as photon-axion mixing [7J[8l[2j[5]. 

Combining direct meauserments of cosmic expansion [30] and SN data [TjJ , and 
assuming an underlying flat ACDM model, we have placed strong constraints on such 
deviations, in particular improving the existing constraint [21] by nearly an order of 
magnitude. 

While allowing for such deviations from transparency, we have revisited the 
reported tension [25j between distance measures obtained from SN and BAO within the 
ACDM model. We have found no discrepancy among Supernova data -even allowing 
for deviations from transparency- , H{z) determinations and high-redshift BAO, 
indicating that the source of the discrepancy may likely be found among systematic 
effects in the modelling of low redshift BAO data (or be attributed to a 2-a statistical 
fluctuation), rather than in exotic physics. 

Finally, since Supernovae data and H (z) data appear to be highly consistent with 
each other, we have combined them and attempted to draw some model-independent 
conclusions about the recent accelerated expansion. We determine the transition 
redshift between deceleration and recent acceleration to be z acc = 0.35^° 13 < 

To conclude, we remark that the comparison among distance measurements 
along the lines of the analysis presented here can also be used to place bounds 
on possible departures from homogeneity and isotropy. There has been significant 
activity recently in trying to use comparisons among distance measurements as a 
powerful tool for probing fundamental assumptions behind cosmological models. For 
example, Clarkson et al [6] suggest to test the Copernican Principle by examining 
the consistency condition: C(z) = 1 + H 2 (DD" - D' 2 ) + HH'DD' = 0, where 
D = (1 + z)~ 1 d,L and a prime stands for differentiation with respect to redshift z. 
This condition arises by solving for the curvature energy density parameter flk from 
the definition of the luminosity distance in a FLRW universe, and then demanding 
that its derivative vanishes. In general, deviations form the Copernican Principle may 
appear as inconsistency between distance measures. The condition above relies on 
the same consistency relation we used here, that is on the consistency relation that 
exists within the homogeneous and isotropic FLRW models between luminosity or 
area distance and the Hubble rate, both as a function of redshift. In particular the 
Clarkson et al condition is not satisfied for radially inhomogeneous models in general. 

While our data are too noisy to perform the test in the form stated above (note 
that C(z) involves first and second derivatives of the data), here, in the language 
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of [B], we have effectively looked at radial deviations from the Copernican Principle 
by testing consistency between different distance measurements {dr,{z) and H(z)). 
We have found no indication of inconsistency among the different datasets/distance 
measures. In general, radial inhomogeneities would lead to anisotropic clustering, 
so, for example, the BAO feature in the two-point correlation function would not be 
isotropic [6J, yielding inconsistency between the inferred d,A and H{z). A similar effect 
is expected in a large class of Bianchi models. We expect that future data will impose 
interesting constraints on these models. 
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